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Abstract 

The Bose-Einstein correlations of photons emitted from a longitudinally expanding 
system of excited matter produced in ultrarelativistic heavy ion collision are studied. 
Two effects found in recent calculations - that the correlation function in longitudinal 
direction exhibits oscillations, and that it takes values below unity - are demonstrated 
to be numerical artefacts and/or the results of inappropriate approximations. Thus, 
the general quantum statistical bounds for the two particle correlation function of a 
chaotic source with Gaussian fluctuations are confirmed. Two different expressions 
for the two-photon inclusive distribution are considered. Depending on which of the 
two expressions is used in the calculation, the width of the correlation function may 
vary by as much as 30%. 
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The study of high energy photons can yield information about the early stages of the 
evolution of hot and dense matter created in ultrarelativistic nuclear collisions. This is of 
special interest in the experimental search for a quark-gluon-plasma (QGP), a new phase 
of matter expected to exist for a brief period of a few fm/c before the majority of the final 
state hadrons are emitted. In particular, information on the space-time extension of the 
excited matter which can be extracted from the two-photon correlation function may be 
used to determine the energy densities reached in the collision. 

Under the title "Photon interferometry of quark-gluon dynamics" , a calculation of the 
Bose-Einstein correlation (BEC) function of photons emitted from a system of longitudi- 
nally expanding hot and dense matter at RHIC and LHC energies was presented in Ref. 
[jij. The results were later extended to include the effects of transverse flow, i.e., to the 
case of a full three-dimenssional hydrodynamic expansion [Q, |j. It was argued Q that the 
BEC's of photons of high transverse momenta (on the order of a few GeV/fm) can be 
sensitive to the presence of a mixed phase and may thus provide a signature for the QGP. 

It is the purpose of this note to point out that two rather surprising features of the 
two-photon correlation function presented in Ref. Jl[ do not reflect genuine properties of 
the source. Rather, they are artefacts of inappropriate approximations used to evaluate 
space-time integrals in Q. To be specific, in Ref. [|J it was found that the BEC function 
in longitudinal direction (a) exhibits oscillations and (b) takes values below unity. As was 
discussed in Ref. [f| , property (b) is inconsistent with the quantum statistical bounds for 
a purely chaotic source. Below, it will be demonstrated that both properties (a) and (b) 
disappear if the integrals are calculated more exactly. 

The BEC function of two identical particles can be written as 

C 2 (ki,k 2 ) = — , (1) 

Pl(fcl)Pl(fc2) 

where fcj {i = 1,2) are the three-momenta of the particles, and P2(ki,k2) and Pi(k) are 
the one- and two-particle inclusive spectra. For the general case of a Gaussian density 
matrix and a purely chaotic source, the current formalism allows to relate all m-particle 
inclusive distributions to a source function w(x, k) which describes the mean number of 
particles of four-momentum k emitted from a source element centered at the space-time 
point x (cf. H and refs. therein). For the one- and two-particle spectra, one has 



Px{k) = / d 4 x w(x,k), (2) 



and 



P 2 {k x M) = Pi(h)P 2 (k 2 ) + J d A Xl d A x 2 w (x l5 h±tl^ w ^ h±tl^ cos(A^A^). 

(3) 
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The expressions (|2|) and (||), which are usually applied in calculations of the BEC function, 
have also been derived in a Wigner function approach in ||. In Ref. Q, however, a different 
form was used for the two-particle inclusive distribution, namely, 

P 2 (k 1 ,k 2 ) = Pi(h)P 2 (k 2 ) + Jd 4 Xl d 4 x 2 w{x 1: ki) w{x 2 ,k 2 ) cos(Ax^AA; M ). (4) 

The form fl§) 

was also derived in Refs. @ as an approximation to the complete, non- 
analytic result. As was shown in 0], the expression (Q) has the disadvantage that in 
certain cases it can lead to results for BEC function which violate the quantum statistical 
bounds for a purely chaotic source, 

1 < C 2 (h,k 2 ) < 2, (5) 

whereas (||) automatically satisfies Eq. (||). 

We have evaluated the above expressions (|3|) and (|j) for the case of thermal photons 
emitted from a longitudinally expanding system of excited matter created in Pb + Pb 
collisions at = 200 AGeV . We use the source function j7|] that was adopted in |IJ, 

, , 9 /2.9k u u' 1 (x) \ ( k u u' I (x)\ 
w(x,k) = const. T{xf ln(— ^-Li + lj exp (-^A 2 ) (6) 

where g is the QCD coupling constant, T{x) is the local temperature and u^(x) the local 
flow velocity. To allow for comparison of our results with those of Ref. jlj], we used the 
same model for the space-time evolution of the dense matter as in O . That is to say, we 
applied Bjorken hydrodynamics || and adopted a massless pion gas equation of state and 
a bag model equation of state for the hadronic phase and the QGP phase, respectively. 
The results presented below correspond to an initial proper time r« = 0.124 fm/c, an 
initial temperature Tj = 532 MeV and a freeze-out temperature Tf = 140 MeV. 

Fig. 1 shows the BEC function]] in longitudinal direction, i.e., for a configuration of 
equal transverse momenta, k%± = k 2 ±. The correlation functions are plotted against the 
rapidity difference Ay = y\ — y 2 . The results presented in Fig. 1 refer to photons emitted 
from the QGP phase only, at a transverse momentum k±_ = 5 GeV/c. 

The curves labeled 1 and 2 were obtained by performing all space-time integrations 

numerically. Curve 1 corresponds to the expression (J3j) and curve 2 to the expression 

(|j) for P 2 (k\,k 2 ). Neither of the two results shows oscillatory behaviour, and they both 

respect the quantum statistical bounds (||). It is, however, interesting to note that the 

width of the correlation function increases by about 30% if (Q) rather than (|3|) is used 

Here and in the following, we consider pairs of photons of equal polarization. Averaging over polariza- 
tions leads to modifications which manifest themselves, among other things, in a decrease of the intercept 
of the correlation function. 



2 



in the calculation. Let us now consider the approximation adopted in Ref. Q. The 
space-time integrals in (^]) — (Q) involve integrations over the space-time rapidities r\\ and 
7]2- In [U, the latter were performed analytically by using a Gaussian approximation of 
the integrand as a function of rji and r] 2 in Eq. (Q). This procedure leads to the results 
displayed as curve 3 in Fig. 1. The corresponding correlation function oscillates as a 
function of Ay and even takes values below one, i.e., it violates the quantum statistical 
bounds (H). A comparison with the "exact" results (curve 2) shows that the Gaussian 
approximation is inappropriate in this case. For completeness, we have also included the 
cosine- approximation proposed in jlj] for small Ay (curve 4 in Fig. 1), 

C 2 (h,k 2 ) » 1 + cos[4k ± Ti sinh 2 (Ay/2)]. (7) 

In Fig. 2 the BEC in longitudinal direction is shown for the sum of all three con- 
tributions (hadron gas, mixed phase and QGP). As was already observed for the QGP 
component alone, at k± = 5 GeV the curves obtained by using the expressions 

d) and (|) 

differ in width by a factor of ~ 1.3. On the other hand, the corresponding two curves for 

k± = 2 GeV are almost indistinguishable. It is noteworthy that while the 5 GeV curves 

are of of approximately Gaussian shape, the 2 GeV curves are distorted and suggest a two 

(or more) component structure. Fig. 3 compares the plasma contributions to the sum of 

the contributions from all phases. For 5 GeV photons the plasma contribution dominates 

whereas for transverse momenta of 2 GeV the shapes of C 2 for the plasma component 

and the sum of all components differ strongly. Indeed, as can seen in Fig. 4 the transition 

from a two-component shape to a Gaussian shape occurs in a narrow range of transverse 

momenta, 2 GeV < k± < 3 GeV. This behaviour can be understood in terms of the 

relative importance of the contributions from the QGP, mixed and hadronic phase. Fig. 5 

shows the production rate of photons at y = as a function of longitudinal proper time r 

for five different values of k±. The two points in proper time where the slope changes dis- 

continuously correspond to the transition from the QGP to the mixed phase and from the 

mixed phase to the hadron gas, repectively. For high transverse momentum photons (e.g., 

k± = 4 GeV), the QGP contribution dominates, and consequently, there appears a single 

Gaussian in the correlation function (cf. Fig. 3). On the other hand, for smaller transverse 

momenta (k± ~ 1 GeV) the contribution from the mixed phase becomes comparable to 

or even exceeds the plasma contribution, resulting in a two-component structure of C 2 . 

Thus, if such a change of shape as a function of k± will be observed it may be considered 

as a piece of evidence for a first order phase transition .0 

2 Of course, other effects can also give rise to a two-component structure. One important example is 
quantum statistical partial coherence |J . 
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In pi, it was argued that a two-component form of the two-photon correlation function 
in transverse direction may signify the presence of a mixed phase and hence, of a first 
order phase transition. In Fig. 6 the BEC function is plotted against the component 
q ou t of the transverse momentum difference (i.e., the component parallel to the transverse 
momentum of the pair), for various values of k±±. The results displayed in the figure 
which were obtained by using the expression (|3|) roughly agree with those obtained in 
Rcf. Q by using (||). As in the case of the correlation function in longitudinal direction, 
one observes the deviation from the Gaussian shape as k±_ decreases, which signifies the 
increasing importance of the contribution from the mixed phase. Note that in contrast to 
the BEC function in Ay, the correlation function in q ou t does show oscillations. However, 
it does not take values below unity, i.e., it does not violate the quantum statistical bounds 

(§)■ 

To summarize, we have demonstrated that two rather surprising properties of the two 
photon correlation function presented in a recent publication are artefacts of inappropri- 
ate approximations in the evaluation of space-time integrals. In it was found that the 
BEC function in longitudinal direction (a) oscillates and (b) takes values below unity. As 
property (b) is inconsistent with general statistical bounds, it is important to clarify the 
origin of this discrepancy. We have shown that both properties (a) and (b) disappear if 
the space-time integrations are performed numerically (i.e., without adopting any analytic 
approximation). On the other hand, we have confirmed that the correlation function in 
transverse direction does exhibit oscillatory behaviour in the out-component of the mo- 
mentum difference, as was observed in Ref. Q. We have considered two different forms 
for the two-particle inclusive distribution and found that the widths of the resultant cor- 
relation functions may differ by up to 30% depending on which one of the two expressions 
is used in the calculation. The change of the BEC function in Ay from a Gaussian to 
a two-component shape with decreasing transverse momentum of the pair may serve as 
evidence for the presence of a mixed phase, and hence, of a first order phase transition. 

Note that the results discussed above refer to a purely longitudinal expansion. While 
there are indications JO]] that at high collision energies the expansion indeed becomes ef- 
fectively one-dimensional, it will be interesting to extend the present considerations to 
the case of a fully three-dimensional hydrodynamic flow. However, this is a task which 
exceeds the scope this brief report. 

This work was supported by the Federal Minister of Research and Technology under 
contract 06MR731, the Deutsche Forschungsgemeinschaft and the Gesellschaft fur Schwer- 
ionenforschung. 
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Figure Captions 



Fig. 1 Two-photon Bose-Einstein correlation C2 as a function of the rapidity difference 
Ay, for photons of transverse momentum k± = 5 GeV emitted in the QGP. The 
curves correspond to different expressions and approximations for the two-photon 
inclusive distribution. 

Fig. 2 Comparison of C2(Ay) calculated from the two expressions (||) and (|j) (see text). 
The results include contributions from all three components (QGP, mixed phase and 
hadron gas). 

Fig. 3 Comparison of the complete correlation function C2(Ay) and the separate contri- 
bution from the QGP phase. 

Fig. 4 Dependence of the BEC function C2(Ay) on the transverse momentum k± of the 
pair. 

Fig. 5 Dependence of the production rate of thermal photons on the longitudinal proper 
time, for different values of the transverse momentum. 

Fig. 6 Two-photon BEC function in transverse direction, as a function of the momentum 
difference in oui-direction, for different values of the transverse momentum k\± of 
one of the two photons. 
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